

! testy MMM

program t1

use MQuick

integer :: i,n
real :: data(10000),d(10000),mean,median,mode,sigma,skew
double precision :: s

i = 0
do
  i = i + 1 
  read(*,*,end=20) data(i)
enddo
20 continue
d = data

call QuickSort(data)
n = size(data)
write(*,*) n
call MMM(n,data, 100.0,mean,median,mode,sigma,skew)

write(*,*) 'n=',n,' mean=',mean,' median=',median,' mode=',mode
write(*,*) 'sig=',sigma,' skew=',skew

n = 10000
s = 0.0
do i = 1,n
  s = s + d(i)
enddo
write(*,*) s/n

end
